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We consider the ground-state properties of the s = i Ising chain in a transverse field 
which varies regularly along the chain having a period of alternation 2. Such a model, 
similarly to its uniform counterpart, exhibits quantum phase transitions. However, the 
number and the position of the quantum phase transition points depend on the strength 
of transverse field modulation. The behaviour in the vicinity of the critical field in most 
cases remains the same as for the uniform chain (i.e. belongs to the square-lattice Ising 
model universality class). However, a new critical behaviour may also arise. We report the 
results for critical exponents obtained partially analytically and partially numerically for 
very long chains consisting of a few thousand sites. 
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The one-dimensional s = ^ Ising model in a transverse field is sometimes called 
a "drosophila" of the quantum phase transition theory since it provides the back- 
ground for the analysis (very often rigorous analytical) of some important aspects 
of quantum phase transitions which occur in more realistic models of real systems 
Q 121 01 ' In what follows we deal with a regularly alternating s = 5 transverse Ising 
chain the ground-state and the low-temperature properties of which have been dis- 
cussed recently in Refs. To be specific we consider N ^ 00 spins ^ governed 
by the Hamiltonian 

i?=-2^«+i+^(f2-(-l)"AI2)<, (1) 

n n 

i.e. we consider the s = Ising chain with the nearest neighbour interaction J = — 2 
in a regularly alternating transverse field of period 2 the two values of which, Sl+Af] 
and ]7 — AQ, are coming into play alternatively. The quantum phase transition 
points for such a chain follow from the condition 6 (J?* + Ai7) (J?* — Ai7) = ±1, 
which immediately yields f2* = ±Vl + Af2^ if AQ < 1, Q* ^ {-\/2, 0, %/2} if 
Af2 — 1 and f2* — ±V Ai?^ ± 1 if Ail > 1. Hereafter we focus on the behaviour of 
different ground-state quantities of the spin chain with Ai7 = 1 in the vicinity 
of critical points H* — \/2 and i7* = since for other values of Ai7 (or for Ai7 = 1 
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and n* = — \/2) the critical behaviour is the same as the one in the vicinity of 

n* = V2. 

A critical behaviour is characterised by a set of exponents which determine pe- 
culiarities of different ground-state quantities in the vicinity of the critical field. 
Thus, the order parameter (i.e. the longitudinal magnetization = w^ni^n)) 
decays as ~ — f^)^ , the equal-time two-point correlation function of x spin 
component behaves as {SnSn+2r) ~ {■^n){^n+2r) j2pjrr exp ^ith the cor- 

relation length ^ ^ The transverse static susceptibility diverges as 

~ iSl—n^l^". Finally, the energy gap vanishes as A ~ ~ | J7 — J?"^ | ''^ , 
where z is the dynamic critical exponent. The quantum phase transition in d space 
dimensions (in our case d = 1) corresponds to the finite-temperature ("classical") 
phase transition 'm d+ z space dimensions. Moreover, from the classical phase tran- 
sition theory we know a number of relations which connect the critical exponents: 
2 - a = 2/3 + -f, 2 - a = P {S + 1), 2 - a = dv, J = {2 - ri) V. 

In the uniform case (A^? = 0) the critical exponents are well known: (3 — ^, 
v = 1, T] = J, a = 0, z = 1. This implies that the quantum phase transition 
belongs to the square-lattice Ising model universality class. In the considered case of 
regularly alternating chain |^ with AJ? = 1 we shall see that the critical behaviour 
in the vicinity of f2* — \/2 remains in the same universality class. However, the 
critical behaviour in the vicinity of i7* = is characterised by a different set of 
exponents. 

We start with the exponents vz and a which can be obtained analytically using 
the Jordan- Wigner fermionization and continued fractions. From Refs. 01 |S] we 
know the elementary excitation energy distribution which permits us to get for the 
energy gap in the vicinity of the critical points, A ^ \i7 — a/2| and A ~ This 
implies vz = 1 and vz — 2 for the critical points H* = \/2 and Q* = 0, respec- 
tively. Moreover, from Refs. 01 E] we know that the ground-state energy contains 
a nonanalytical contribution, eo ~ (i7 — \/2) ln|i7 — \/2| and eg ^ f2'*ln|J7|, in 
the vicinity of ]7 = V2 and = 0, respectively. After differentiating we conclude 
that = %^ exhibits a logarithmic divergency at = \/2 and remains finite 
at i7* = where, nevertheless, exhibits a logarithmic divergency. Therefore, 
a = at 12* = \/2 and a = -2 at /?* = 0. 

We go on to the exponents /3, v and rj which can be calculated numerically 0. 
Knowing the spin correlation function {SnSn+2r) '^^^ obtain the on-site magne- 
tizations mfj — y/limr^oo{SnSn+2r) ' n = 1,2 and = i {mf -|-to|). Assuming 
for {SnSn+2r) ^ the long-distance behaviour ~ (2r)~''exp as r ^ oo 

we can find the correlation length ^ and the exponent rj. To determine (Fig- 
we take TV = 2000, 4000, 5400, n = f , 2r = f . To determine ^ and r] (Figs. (HO 
we take N ^ 2000 and n = 920, N = 4000 and n = 1920, N = 5400 and n = 2620, 
and 2r = 162, 170, 178, 186, 194. 

Let us comment on the the obtained results. From the data reported in Figs, 
n 12 El we may conclude that the critical behaviour in the vicinity of il* — \/2 
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Fig. 1. The order parameter behaviour in the vicinity of critical points (A'^ = 2000: 
squares, N = 4000: diamonds, N = 5400: circles). The solid lines in panels a, b and c 

represent the dependences m"" = 0.5 (-s/S ± J?)^ and m'^ = 0.5|i?|3 (a), m"^ = 0.5Oi (b) 

and = 0.5 (\/2 - r?) ^ (c). 
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Fig. 2. The inverse correlation length behaviour in the vicinity of critical points (A'' — 2000: 
squares, N = 4000: diamonds, N = 5400: circles; diamonds coincide with circles in the 
scale of panels b and c). The solid lines represent the dependences | = 3\\/2 — if 
n <V2 {a and c), | = V2\V2 - i7| if J? > \/2 (a and c) and | = (a and b). 



is characterised by the set f3 = 0.125, u = 1 and rj = 0.25. Contrary, in the 
vicinity of = the corresponding exponents are f3 = 0.25, ly = 2 and 77 = 0.25. 
Obviously, we expect such N ^ 00 values of /3, v and rj from the finite-A^ data 
fohowing the tendencies in their change as N increases. Combining these findings 
with the analytical results for vz and a we get two sets of critical exponents. 
Namely, P = 0.125, v = 1, t] = 0.25, a = 0, vz = 1 in the vicinity of f^* = \/2 and 
/3 = 0.25, v = 2,r\ = 0.25, a = -2,vz = 2 in the vicinity of 17* = 0. The former set 
corresponds to the square-lattice Ising model universality class whereas the latter 
one corresponds to a new universality class. 

To summarize, we have studied the quantum phase transitions in the s = | Ising 
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Fig. 3. The power-law exponent r] in the vicinity of critical points (A'' = 2000: squares, 
N = 4000: diamonds, N = 5400: circles (they cover diamonds in the scale of panels b and 
c)). T] tends to 0.25 at the critical points very slowly with increasing N . 



chain in a regularly alternating transverse field of period 2, i7± AJ7. We have found 
that for a small strength of regular nonuniformity (controlled by Al?) the ground- 
state phase transition between the Lsing phase (for small fi) and the paramagnetic 
phase (for large Q) is modified only quantitatively, i.e. the value of the (two) critical 
fields increases. The critical behaviour remains of the square-lattice Ising model 
universality class. If Af? is equal to a half of the exchange Ising interaction, an 
additional quantum phase transition point appears at J7 = 0. We have found a 
strong evidence that this critical point is characterised by the new set of exponents. 
The other (two) critical points which correspond to the transition between the Ising 
and paramagnetic phases are of the square-lattice Ising model universality class. 
For larger Af? the chain exhibits four quantum phase transition points between 
the Ising and paramagnetic phases which belong to the square-lattice Ising model 
universality class. We should stress that a part of the exponents have been obtained 
numerically. Although, in principle, the spin correlation functions required for the 
estimation of these exponents can be calculated analytically as in the imiform case, 
to our best knowledge, it has not yet been done. We have left this problem for 
future studies. 

The authors are thankful to Johannes Richter and Oles' Zaburannyi for collaboration 
in the earlier stage of this study. O. D. has enjoyed the hospitality of the guest program 
of the Max-Planck-Institut fiir Physik komplexer Systeme, Dresden in the spring of 2004. 
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